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TECHNICAL MEMORANDUM 


HEAT TRANSFER AND PRESSURE DROP IN CONCENTRIC ANNULAR FLOWS 

OF BINARY INERT GAS MIXTURES 

1. GENERAL PROJECT OVERVIEW 


NASA Marshall Space Flight Center (MSFC) builds components for nonnuclear thermal- 
hydraulic testing of fast spectrum reactor cooling systems for advanced, high-powered, deep-space 
probes and surface power systems. One option under consideration is to couple a gas-cooled reactor with 
annular flow passages to closed-loop gas turbine engines. Such engines use mixtures of heavy and light 
inert gases (typically helium (He)-xenon (Xe)) to reduce core and heat exchanger size and to decrease 
the number of compression stages. 1 A notable feature of mixtures of gases with differing molecular 
weight is that they have unusually low Prandtl numbers (Pr). Figure 1 compares the Prandtl number 
of three different gas mixtures as a function of the lighter gas fraction (in all cases, He). This activ- 
ity focuses on an experimental investigation of Prandtl number effects using test sections with various 
annulus ratios. 



Figure 1. Prandtl number versus mole fraction for three binary 
gases near room temperature. 


1.1 Objective 


The principle objective of this test program is to investigate heat transfer coefficients in 
centrally heated annular ducts at low Prandtl numbers. The immediate need is to obtain data that will 
be used for comparison to analytical models and existing heat transfer correlations to guide reactor 
design. Specifically, the data will address two issues previously identified as the low Prandtl number 
heat transfer effect and the large channel decrement effect. These effects could potentially impact reactor 
design. Specific supporting objectives are as follows: 

• Provide Nusselt number (Nu) data using He-argon (Ar) gas mixtures for various annular duct 
geometries over a range of operating conditions. 

• Compare results to existing correlations and recommend correlation for design use. 

The completion of this test phase will lead to tests to investigate heat transfer at prototypic 
temperatures and pressures using He-Xe gas. 

1.2 Background 

Current reactor design work supporting Project Prometheus is focused on a gas-cooled reac- 
tor using an He-Xe coolant mixture. Prototypic thermal hydraulic data are required for ongoing design 
development. These data include pressure drop and heat transfer characteristics for the geometries and 
system conditions currently under consideration. In addition, two issues have arisen that are affecting the 
analytical design of the reactor: 

• Lack of heat transfer correlations for gas mixtures with Prandtl numbers of ~0.2. 

• Potential correlation inaccuracy in evaluation of coolant passages with large thickness/diameter. 

The range of gas molecular weight under consideration is 20 to 40 g/mole. The baseline He-Xe 
core inlet temperature is 890 K and the outlet temperature is 1,150 K. The core inlet temperature could 
vary between 810 and 900 K. The core outlet temperature could vary between 1,050 and 1,150 K. The 
baseline pressure for the reactor design is 2 MPa, with a range of 1.38 to 2.5 MPa pressure being consid- 
ered. In support of other activities, MSFC has assembled a single-channel test rig for characterization 
of pressure drop and heat transfer rates through an electrically heated gas channel.- Nitrogen (N 2 ) 
flowed through this rig in an open, single-pass arrangement. MSFC has been tasked to modify the 
single-channel test article for heat transfer and pressure drop characterization of various gas mixtures. 
Test sections of immediate interest include annular pin, circular ducts, and axial flow open lattice con- 
figurations. Test data generated will include pressure drop, and surface and bulk temperatures sufficient 
to derive the heat transfer coefficient and friction coefficient as a function of system conditions and local 
position. First priority will be given to the existing annular pin test configuration with a smooth fuel tube 
insert. All annular tests will have a heated core tube and adiabatic outer tube. Tests may also be con- 
ducted on an electrically heated tube to determine correlations for in-tube (nonannular) flow. It follows 
then that the annulus ratios under consideration may range 0<r*< 1. 
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1.3 Significance of Results 


The film temperature drop is a large contributor to the overall fuel system temperature. Heat 
transfer test results were intended to guide the reactor methods development in support of gas-cooled 
reactor design. The results to be obtained have the potential to dramatically change the reactor thermal 
design (30%-50% heat transfer decrement). Recent computational fluid dynamics (CFD) analyses have 
indicated that increasing coolant passage thickness can result in a lower heat transfer coefficient than is 
predicted by correlations. Testing of various sized coolant passages will investigate this analytical trend. 
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2. LITERATURE REVIEW 


2.1 Thermophysical Properties 

Problems involving fluid flow and heat transfer of gas mixtures requires knowledge of the 
viscosities, thermal conductivities, and Prandtl numbers as a function of composition and temperature. 
Direct measurements of such properties for pure species have been, in most cases, conducted over a 
wide range of temperature for rarefied gases. Transport properties for pure gases, such as viscosity and 
thermal conductivity, can be found using polynomial curve fits as a function of temperature. Higher 
order polynomial curve fits for pure gases have been made from accepted values found in Touloukian. 3 ’ 4 
These properties are found in appendix A. 

The Prandtl number is defined as Pr=[ic p /k. For a monatomic ideal gas, the gas constant is 
R = c p -c v and the ratio of specific heat, y=c /c v , is 5/3, so the specific heat capacity at constant pressure 
is c | =R/M\( 1-3/5)“ 1 . For a mixture of ideal monatomic gases, the heat capacity is c p M = m ] c p j 
+ (l-m 1 )c 2 ’ where the mass fraction is mj = n 1 M 1 /n 1 M 1 + ( 1 ■ 

Properties for mixtures have been measured near room temperature and, in some cases, at 
select temperatures, all for rarified gases. Mixture viscosities and conductivities can be calculated from 
approximations derived from rigorous kinetic theory of gases as applied to a realistic intermolecular 
force laws. 5 The relative thermal conductivity and viscosity of binary gas mixtures cannot be obtained 
by simple averaging of pure gas values. More complex methods of calculation of transport properties of 
binary mixtures of inert gases have been devised. The viscosity of a mixture of gases can be calculated 
by an approach attributed by various authors to Sutherland and found in Vanco. 6 For a binary mixture, 
the viscosity equation can be expressed as 


where 
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The pure gas values are taken to be the values determined by the polynomial curve fits outlined 
in appendix A. The thermal conductivity of a binary mixture of monatomic gases is calculated with 
a method using a similar Sutherland formulation: 


where 


and 
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Figures 2 and 3 compare the viscosity and thermal conductivity as a function of He mole fraction 
for an He-Ar mixture at 291.2 K. Such data for He-Xe exists at only a few discrete temperatures. Similar 
plots for He-Xe mixtures are shown in figures 4 and 5. 



Figure 2. Viscosity of He-Ar mixtures from data found in Touloukian 4 (symbols) 
compared to gas mixture viscosity relation in equation (1) (line). 

Gas temperature is 291 K. 
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M (kg/m-s) 



Figure 3. Thermal conductivity of He-Ar mixtures from accepted values 

found in Touloukian 3 (symbols) compared to gas mixture thermal 
conductivity relation in equation (4) (line). Gas temperature is 291 K. 



Figure 4. Viscosity of He-Xe mixtures from 

data found in Touloukian 4 (symbols) 
compared to gas mixture viscosity 
relation from equation (1) (line). 



Mole Fraction of He 


Figure 5. Thermal conductivity of He-Xe 
mixtures from accepted values 
found in Touloukian 3 (symbols) 
compared to gas mixture thermal 
conductivity relation from 
equation (4) (line). 
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2.2 Pressure Drop in Annular Channels 

The present study undertakes measurement of friction factors along smooth annular channels 
as an ancillary to measurement of dimensionless heat transfer. The friction factor,/, for smooth tubes is 
strictly a function of Reynolds number, Re = pVd h p -1 , where for annular cross sections, the length scale 
is defined as four times the cross-sectional area occupied by the fluid divided by the wetted perimeter 
that the fluid exerts skin friction. Fluid friction in annular ducts has been previously characterized in the 
literature. Walker measured fluid friction factors in smooth concentric annuli in the viscous, transition, 
and turbulent range. 7 The modified hydraulic radius concept was found adequate to correlate observa- 
tions. For turbulent flow through a circular tube annuli, Rohsenow recommends use of a Blasius-type 

friction factor formula: / = 0.085RcT°' 25 for 6,000 < Re < 300, 000. 8 

2.3 Heat Transfer in Annular Channels 

Although heat transfer through annular channels has been extensively examined by various 
groups (for instance, for smooth annuli and for ribbed annuli in He), highly accurate benchmarked heat 
transfer correlations in the 0.2</V<0.6 range have not been systematically developed. Taylor, Bauer, 
and McEligot document several correlations that could potentially be used in the evaluation of gases 
with Prandtl numbers of ~0.2. 9 However, few of these correlations were created for use with such gases. 
In addition, the dimensionless heat transfer coefficients produced by these correlations agree poorly with 
one another. Testing of the specific geometries under consideration will help guide the designer in the 
choice of an appropriate correlation. 

Investigators that have examined gas flows through channels of similar geometry or Prandtl num- 
ber condition include Kays and Leung, who considered analytically the heat transfer coefficients of the 
smooth annulus for a wide range of Reynolds number, annulus ratio, and Prandtl number down to 0.6. 10 
Wilson analyzed heat transfer for fully developed turbulent flows in concentric smooth annuli. 11 Taylor, 
Bauer, and McEligot, and later McEligot and Taylor, considered internal-forced convection in tubes of 
low Prandtl number gas mixtures. 12 ’ 13 They found that two Nusselt relations best captured the Prandtl 
number dependency for mixtures in the range 0.2</V<0.7, the Kays correlation: 14 ’ 15 

Nu =0.022 Re 0 8 Pr ' 0 6 , (7) 


and the Petukhov Popov correlation, 16 


(£/8 )RePr 

K^) + K 2 (Pr)(U 8 ) 1 , 2 / 2 / 3 - i )' 


(8) 


where 

£ = (l. 82 log (Re)- 1.64 ) -2 


and 

= 1 + 3.4^ K 2 (Pr) = 11.7 + 1.8/V -1 / 3 . 


(9) 

( 10 ) 
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Figure 6 compares the Kays and Petukhov correlations. Deviation between the two correlations 
is greatest at low Reynolds number. Other correlations developed for liquid metals or used in heat trans- 
fer applications at higher Prandtl number (Colburn and Dittus-Boelter equations 17,18 ) were found defi- 
cient. Later, the same research group examined turbulent Prandtl number in the near wall region for low 
Prandtl number flows. 13 



Figure 6. Comparison of Kays correlation 14,15 (symbols) to Petukhov Popov correlation 16 
(solid lines) as a function of Prandtl number and Reynolds number. 


Heat transfer and pressure drop from augmented surfaces have been also studied. Maubach mea- 
sured pressure drop through annuli and examined the effect of roughness on turbulent velocity profile 
of gas-cooled annuli. 19 Dalle-Donne derived correlations for gas-cooled fuel rods with and without rec- 
tangular ribs. 20 Ichimiya conducted tests on the effects of several roughness elements on insulated walls 
for heat transfer from the opposite smooth heated surface in a parallel plate duct. 21 Takase considered 
turbulent heat transfer performance for He coolant of a fuel rod with spacer ribs for a high-temperature 
gas reactor. 22 

Of primary importance to this investigation is an understanding of the functional dependences 
in the Nusselt correlation for fully developed low Prandtl number flows inside annular ducts. For low 
Prandtl number gas mixtures, the thermal boundary layer develops readily while the momentum bound- 
ary layer development lags, provided there exists no unheated starting length. Nusselt correlations 
typically assume constant thermophysical properties. This assumption may not hold well for the large 
thermal gradients encountered in the Brayton heat exchangers or reactor core. Accurate gas properties 
will be needed as a function of temperature and composition. One approach to separate the first-order 
Nusselt correlation dependencies from fluid property dependencies is to conduct initial tests near 
ambient temperature where gas mixture properties are well established and constant properties can be 
assumed. Once a constant property Nusselt relation at low temperature has been established, tests can 
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be conducted at high temperature and power to separate nonconstant fluid property effects. Later, tests 
at near prototypic reactor conditions can capture the effects of nondeveloped flow conditions associated 
with nonuniform axial heating profiles. 
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3. TEST CONDITIONS 


3.1 Test Summary 

Table 1 shows nominal dimensions for the channels under evaluation for the single-channel test. 
Gas in the annulus is heated along the inner tube wall with a uniform heat flux. For the first series of 
tests, one inner tube size is used and the jacket diameter is changed to produce different annulus ratios. 
The length of the annulus and heated region of each test section is ~18 in. The inlet conditions are con- 
trolled so that thermal- and hydrodynamic-developed flow exists along a significant fraction of the test 
section length. The first phase of testing will use an inner tube size of 0.625 in outside diameter (OD). 
The jacket inside diameter (ID) will be varied (0.750, 0.815, 0.875, and 1.250 in inside diameters) to 
obtain the different annular ratios. The annular flow channel will be finished to <0.8-pm rms to ensure 
hydrodyamically smooth surfaces. 


Table 1. Channel dimensions for single-channel test. 


Test Section 
Designation 

Geometry 

Type 

Inner Tube 
OD (in) 

Jacket 
ID (in) 

Annulus 

Ratio 

AF1 

Annulus 

0.625 

0.750 

0.833 

AF2 

Annulus 

0.625 

0.815 

0.767 

AF3 

Annulus 

0.625 

0.875 

0.714 

AF4 

Annulus 

0.625 

1.250 

0.500 


Figure 7 shows traces of the four annular duct geometries. Gas in the annulus is heated along 
the inner tube wall with a uniform heat flux. Care must be taken in the positioning of the inner tube rela- 
tive to the outer jacket tube to minimize eccentricity of the flow space, as it has been shown by Judd and 
Wade to have a significant effect on local Nusselt numbers. 8 As will be discussed later, proper position- 
ing is especially important for the test section labeled AF1. 



Figure 7. Cross-sectional traces of four annular flow passages. 
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3.2 Test Matrix 


Testing of the single channel will begin with a pure gas, followed by tests with an He-Ar mixture 
that produces Pr ~0.4. Initial tests will be near room temperature and focus on heat transfer measure- 
ments for constant property flows. Variables in the test conditions include gas composition (Prandtl num- 
ber) and Reynolds number. Data will be taken near test section thermal equilibrium. Inlet temperature 
to the single channel will be near room temperature. Table 2 gives the range of conditions for the initial 
He-Ar test series of each test section. 


Table 2. Nominal tests. 


Test 

Pr 

Re 

1 

0.67 

Transition 

2 

0.67 

7,000 

3 

0.67 

15,000 

4 

0.67 

22,000 

5 

0.67 

30,000 

6 

0.40 

Transition 

7 

0.40 

7,000 

8 

0.40 

15,000 

9 

0.40 

22,000 

10 

0.40 

30,000 


Variables in the test conditions include He-Ar mixture (Prandtl number) and Reynolds number. 
Data are to be taken at steady- state conditions for uniform heating of each test section. Inlet temperature 
to the single channel will be initially near room temperature. Mass flow rate and power is adjusted to 
create a specified Reynolds number and temperature rise across the test section. For each gas mixture 
under investigation, the first step will be to identify the transition point from laminar to turbulent flow. 

Required data include total power input to test section, mass flow rate, heated surface tempera- 
ture as a function of axial position, bulk coolant temperature as a function of axial position, pressure 
at inlet and outlet, incremental pressure drop along the heated region, and calculated Nusselt number 
for each test condition. 

Initial testing will focus on test section AF2 for pure gas testing with Ar followed by He-Ar test- 
ing. As a benchmark, pure gas test results will be compared to exiting data in the literature for pure gas 
heat transfer. After successful acquisition of Prandtl number-sensitivity data with AF2 are complete. 
He-Ar testing of the remaining test sections will be completed. To ensure data accuracy, repeat data 
points will be taken for a subset of the test matrix several days or longer after the initial data acquisition. 
Also, single-sample uncertainty calculations will be made to identify significant error sources. Heat 
loss testing will be performed to quantify system losses. Zero-power pressure drop testing will be per- 
formed. Energy balance calculations for the system will also be performed to verify consistency in the 
data. Additional data quality information will be produced, including verification of the concentricity 
of the inner tube and the outer jacket, methodology of ensuring fully developed flow (with supporting 
evidence), and a procedure for identification of the transition Reynolds number. The expected instrumen- 
tation for the first phase of test includes temperature, pressure, pressure drop, gas composition, power 
to test section, and gas flow rate. 
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The second phase of this test will require a closed He-Xe loop operated at prototypical condi- 
tions— pressure, temperature, flow, and gas composition. Geometry investigations will be focused on 
the geometry selected at the end of the preconceptual design phase. (Alternate geometries under consid- 
eration include circular duct and open lattice flow.) 


12 


4. EXPERIMENTAL APPARATUS 


The single-channel tests make use of facilities developed for a heat pipe heat exchanger compo- 
nent test. 2 A number of modifications will be made to these facilities to better meet the expanded scope 
of this test. Figure 8 is a photograph of the existing test facility housed in building 4655 at MSFC. To 
accomplish the outline testing, a portion of the floor space within the building 4655 high bay was desig- 
nated to house the setup and operation. An area at the rear of the laboratory with an approximate rectan- 
gular footprint of 15 x 25 ft was set aside. Geometric constraints were defined for each of the primary test 
support system so a revised laboratory layout could be generated. These systems include the following: 

• Test section. 

• Test loop. 

• Inert gas source. 

• Flouse cooling water loop. 

• Power distribution system. 

• Instrumentation and control. 



Figure 8. Existing building 4655 test facility. 
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The location of building 4655 provides close access to both facility water and power interfaces 
while minimizing interference with normal day-to-day operation. The original facility was configured 
for single-channel tests with pure N 9 in a single -pass flow arrangement. While an open flow arrangement 
is appropriate where the gas supply is plentiful, a closed loop better conserves more expensive gases 
such as Ar or Xe. 

Figure 8 is a photograph of the floor space in the summer of 2005 just before hardware setup 
began. The existing single-channel rig and the racks to house the data acquisition/control hardware 
and computers are shown on the lower left side. 

4.1 Gas Circulation Loop Design 

Gas is supplied to the annular test section from a closed circulation loop. The loop and test sec- 
tion as now configured can handle inexpensive gases such as He, Ar, or N 2 . Use of more expensive gases 
like Xe will require a leak-free compressor and a cold trap or other method to recover gas from the loop. 
Most of the components for a loop testing low-temperature gas can also be used for tests at high tem- 
perature. The high-temperature design requires components to preheat and cool the gas on either end 
of the test section. 


4.2 Low-Temperature Loop 

Figure 9 shows a schematic of a flow loop for tests near room temperature. The flow loop con- 
sists of a main loop and a nested bypass loop. The low-temperature loop is built from austenitic stainless 
steel and makes use of Conflat flanges or welded joints, whichever is convenient. For the initial series 
of tests with a reciprocating compressor, the flow in the bypass loop may be as much as 10 times that 
in the main loop. The test section can be evacuated with a roughing pump located at the exit of the com- 
pressor. Gas is supplied to the loop from premixed bottles with composition known to be better than 
a percent. A residual gas analyzer is used to confirm composition of the circulating gas. Loop pressure 
is controlled with a Tescom proportional pressure controller that will set the system to a nominal 
100 psia. An in-line pressure relief value is connected near the gas supply and compressor exit. 
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Figure 9. Single-channel, low-temperature test loop. 


On exiting the compressor, compression heat is removed from the gas with a heat exchanger con- 
nected to a house water circuit. This heat exchanger cools the gas to near room temperature. The flow is 
then split between the main loop and the bypass loop. The proportion of gas to each loop is established 
with a variable position valve located in the main line that acts as a throttle and as another valve in the 
bypass line. Gas enters the main loop. In the previous open-loop system, mass flow rates were measured 
with a thermal meter that operates with a calibration that has a complex dependence on gas composition 
and local heat transfer rates. Such a meter would be difficult to use for gas with variable composition. 

For the current setup, gas passes through a Meriam Instruments model 50MH10-2 laminar flow 
element that measures volumetric flow rate. Calibration for the laminar flow meter merely requires 
knowledge of gas density. The gas may optionally enter a Poly Science model 6706P chiller capable 
of bringing the gas to =263 K. The gas then moves through the test section where it is heated along 
the annular channel with power supplied from two high direct-current bus bars. 

For the initial test series, the temperature rise across the test section will be held constant by 
controlling the power to the test section at a given mass flow rate. Gas inlet and exit temperatures are 
measured as well as wall temperature along the inner flow channel and pressure along the outer surface 
of the channel. Heat loss to the environment is controlled with insulation. Gas exits the test section and 
passes through a variable position valve before returning to the bypass loop where it is compressed. 
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Given a suitable compressor or blower, the low -temperature loop can use either He-Ar or He-Xe 
mixtures. The loop lines were sized with 2-in pipe so as not to exclude the possibility of compressing 
the gas with a turbo blower. One drawback of larger lines is the added volume of gas required to charge 
the system. This problem is especially acute for costly Xe. To partly compensate for this, line lengths 
will be kept short throughout and operating pressure kept to the lowest possible level consistent with 
compressor design specifications. 


4.3 High-Temperature Loop 

Figure 10 shows one possible embodiment of a high-temperature gas circulation loop. The 
bypass portion of the circuit is identical to the low-temperature version. However, the main loop has 
a number of changes from the low-temperature design, including movement of the laminar flow meter 
from the test section inlet to downstream of the test section aftercooler. For parts seeing high tempera- 
ture, metal sealed flanges are impractical, so all welded construction will be mandatory in the main loop 
near the test section. Also, use of superalloy material may be required if the test section must operate 
much above 700 °C. On exiting the bypass, the gas enters a preheater that brings the gas to the test sec- 
tion inlet temperature. The gas then enters the test section where it is electrically heated in the annulus. 
A spider and outlet centering ring, made of alumina, electrically isolate the inner tube from the jacket. 
Heat loss to the environment is controlled with insulation and electrically heated guard heaters distrib- 
uted axially along the test section. The gas exits the test section where it is cooled to near room tem- 
perature with a heat exchanger (aftercooler) connected to a house circuit. Once cooled, the gas can pass 
through lines with conventional construction taken from the low -temperature loop. 
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4.4 Compressor Options 


Gas circulation through the loop may be achieved in a number of ways, depending mostly on 
factors unrelated to operating efficiency. Ideally, the compressor would be capable of handling gases 
ranging from pure He to pure Xe at high pressure. An industrial, oil-free, high-pressure gas compressor 
manufactured by RIX Industries, Benicia, CA, is available as part of the existing test facility. This large 
compressor is capable of circulating up to 0.1 kg/s of pure He from 300 to 375 psia. Such capacity far 
exceeds what is required for the single-channel test apparatus. However, the gas leakage rate through 
the RIX compressor’s labyrinth seals is quite large, requiring the change of several K bottles of gas 
during an 8-hr operating period at 100 psia. Use of a compressor that leaks at this rate is an entirely 
practical option for premixed He-Ar that costs ~$200 per K bottle. Premixed K bottles of Xe cost over 
$20,000 each, making the RIX compressor a prohibitively expensive option for an He-Xe system. 

The RIX compressor will be used for the initial tests with He-Ar. Another leak-tight option must be 
devised for the He-Xe test series. Candidate compression systems for use with He-Xe include (1) a Metal 
Bellows model MB-601 pump contained inside a pressurized vessel, (2) a Haskel gas booster pump 
(fig. 11), and (3) a Fuji ring compressor contained inside a pressurized vessel. 



Figure 11. Gas booster pump procured from Haskel for compression of He-Xe mixtures. 
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4.5 Test Section Design 


A CAD rendering of the single-channel test section is shown in figure 12. More detailed draw- 
ings of the test section are shown in appendix B. The test section consists of an inlet and exit manifold, 
jacket tube, and inner tube with associated power feeds, and temperature and pressure sensors. The test 
section is easy to insulate and permits the inclusion of guard heaters to reduce radial thermal gradients 
for a high-temperature design. The low-temperature design can be transformed to a high-temperature 
unit by replacing flanges with welds. Only minor modifications to the test section design are required 
to change the jacket dimensions for tests at different annulus ratios. The modular design of the test sec- 
tion allows multiple vendors to fabricate subcomponents. The test section may be oriented vertically 
or horizontally, permitting the effect of natural convection to be examined. By directly heating the jacket 
tube, and by incorporating an appropriate electrical grounding scheme, this design may be modified 
for circular tube geometry. 



Figure 12. Single-channel test section. 


In its present embodiment, the test section is instrumented with 10 embedded 0.01 -in-diameter 
type K sheathed thermocouples with special limits of error with accuracy on the order of +1.1 °C. These 
thermocouples are electrically insulated from the sheath to avoid ground loops between the data and 
power systems. The sheath of each thermocouple is also insulated with polyolefin shrink tube in places 
where the thermocouple might make contact with the jacket or manifold. Additional sheathed and insu- 
lated thermocouples are placed at three locations at both the inlet and exit of the test section to measure 
the inlet and exit gas temperature. 

Pressure taps are placed at eight locations along the test section as well as at the test section 
entrance and exit. Differential pressure measurements are made along the length of the test section with 
a multichannel pressure scanner. Manufactured by Pressure Systems, the NetScanner™ model 9116 
consists of up to 16 silicon piezoresistive pressure sensors, all referenced to a common port with +0.05% 
full-scale accuracy. For this test, model 9116 is instrumented with 7-kPa, full-scale pressure transduc- 
ers and is capable of operating at line pressures up to 4 MPa. Eight pressure taps are located along the 
length of the flow annulus with the fourth tap from the test section entrance serving as the reference port. 
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4.6 Mechanical Alignment 


For well-controlled heat transfer measurements, the eccentricity of the inner tube with respect 
to the jacket tube is a key variable to control, especially when the annulus ratio exceeds ~0.8. Annulus 
eccentricity is the distance between the centerlines of the inner and jacket tubes divided by the mean gap 
between the tubes. Eccentricity can affect local Nusselt numbers as shown in Judd and Wade. 8 Eccen- 
tricity is estimated based on measurements of total indicated runout of the inner and jacket tubes for the 
AF2 test section, which are 0.0085 and 0.005 in, respectively. Using this estimate, the eccentricity of the 
assembly for test section AF2 is approximately 

2e * 2 (fwj-) +K) 2 ( (0.0085 ) 2 + (0.005 ) 2 j 

Sc — ~ — — — 

Di-do Dj - d a 0.815-0.625 ( n ) 

e = 0± 0.104 (20 to 1). (12) 

Figure 13 shows the relation between eccentricity and local Nusselt number at various annulus 
ratios. For annulus ratios <0.8 and eccentricity <0.1, an approximate linear relation between Nusselt 
number and eccentricity can be used: Nu/Nu conc ~ 1 -e/2. Outside these ranges, the relation becomes 
nonlinear and an evaluation must be made for each annulus ratio. For an annulus ratio of 0.83, the poly- 
nomial fit for the eccentricity factor is Nu/Nu conc ~ - 1 00 U + 5e 2 -e+l. 



e 


Figure 13. Influence of eccentricity and annulus ratio on turbulent-forced 
convection heat transfer in annular passages. 


Centering rings produce concentric alignment between the inner tube and the jacket. A spider 
located at the inlet of the test section maintains alignment once the centering ring at the test section inlet 
is removed. The outlet centering ring is retained after test section assembly, holding the inner tube at the 
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exit of the annular region. Both the spider and the outlet centering ring are made from G-10 fiberglass 
(for low-temperature design) or alumina (for high-temperature design) to isolate the electrically heated 
inner tube from the jacket body. These centering features collectively permit quick changeout and posi- 
tioning of the inner tube. 


4.7 Inlet Manifold 

A honeycomb core manufactured by Benecor, Inc., Wichita, KS, is located inside the inlet mani- 
fold. This core, shown in figure 14, is a right circular cylinder (4 in diameter by 4 in long) made up of a 
0.002-in-thick sheet that forms cells 0.1 in on a side. A hole is cut into the center axis of this core to pass 
a power lead between the inlet manifold and the inner tube power terminal. This core reduces the Reyn- 
olds number of the gas passing through it, providing a definable boundary condition near the inlet of 
the annular region. The gas then enters a short open region, and then passes into the annular test section 
through a rounded entrance region with a 0.2-in radius. 



Figure 14. Floneycomb structure used in manifold near test section entrance. 


4.8 Test Section Inner Tube 

The test section inner tube forms the core of the annulus; as such, it must simulate the dissipation 
of fission heat. To produce usable heat transfer correlations for fully developed flow conditions, the local 
heating rate should occur at a known and preferably spatially uniform rate. Table 3 gives the test sec- 
tion power requirement for various gases and annulus ratios to maintain Re = 30,000 with a 271 K inlet 
temperature and a 311 K exit temperature. Use of gas mixtures between Fle-40Ar and Fie- 100 Ar span 
0.4 < P/- <0.7. This covers the same Prandtl number range as He-OAr to FIe-39Ar, and owing to the lower 
heat capacity of Ar, requires less power to heat the test section than with He-rich mixtures. 
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Table 3. Test section power requirements (W) for He-Ar mixtures at = 30,000. 


Test Section 

Jacket ID (in) 

He-0%Ar 

He-40%Ar 

He-100%Ar 

AF1 

0.750 

3,283 

848 

371 

AF2 

0.815 

3,438 

888 

389 

AF3 

0.875 

3,582 

925 

405 

AF4 

1.250 

4,477 

1,156 

506 


Two approaches to heating the test section inner tube were considered. One approach uses a 
graphite resistance heater, recently developed at MSFC for other thermal simulator work. A version of 
these graphite heaters has a circular cross section with a diameter of 0.3 in, split along the axis, forming 
two semicircular halves along the heated length. These semicircular halves form a ‘U’ shape, with the 
bottom of the ‘U’ a transition to a circular (nonsplit) cross section at the unheated end opposite the leads. 
Each graphite semicircular half is separated by electrically insulating spacers. The heater ends at the top 
of the ‘U’ connect to power leads. Alumina rings isolate the graphite conductors from the walls of the 
inner tube assembly. In the current embodiment, this heater is inserted into an annular inner tube having 
a 0.625-in OD and a 0.457-in ID. Direct current, applied across the leads located at the test section exit, 
dissipates heat volumetrically. 

Use of graphite heaters for this particular application has a serious drawback. Heat transfer 
between the graphite heater and the inner tube is by radiation, and, in the presence of gas, natural con- 
vection. Proper operation of the test article requires heat to be uniformly distributed on the surface of 
the inner tube. Heat transfer by conduction and radiation are consistent with this requirement. However, 
natural convection of gas between the hot graphite heater and the cooler inner tube potentially induces 
nonuniform distribution of heat on the surface of the inner tube and can have strong dependence on 
the orientation of the test section. Convection can be suppressed by evacuating the volume between 
the heater and the inner tube. Evacuating the inner tube volume leaves only radiation and conduction 
to remove heat from the leads. Unfortunately, as shown in figure 15, the leads under these conditions 
become quite hot, approaching 1,000 K, thus requiring a robust electrically insulating, high-temperature 
hermetic seal. 
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Figure 15. Temperature of a graphite heater cooled by radiation and conduction alone. 

To avoid the vexations associated with high-temperature seals, direct resistance heating of the 
inner tube is now being pursued. In anticipation of this contingency, a 1,400- A, 15-kW power supply 
was purchased from Lambda Americas, Neptune, NJ. Copper bussing attached to each end of the inner 
tube that runs to power terminals are located on opposing manifolds. The sheaths of thermocouples that 
run between the inner tube and the manifolds are electrically isolated from the manifold with polyolefin 
shrink tube. The heating rate is then established, assuming electrical power is uniformly dissipated to the 
inner tube with homogeneous composition and cross section. Additional heating or cooling can occur 
from parasitic losses across insulation that couples the test section gas and the environment: 


An ~2. 3-in-thick layer of insulation, having 0.068-W/m-K thermal conductivity, placed on a 1-in- 
diameter jacket tube is sufficient to keep the second parasitic term <4 W (without other loses). 

Figure 16 shows a method of producing axial holes along the 25-in length of the inner tube. 


0.015 x 0.0 15 in are cut or milled on the outside diameter of a 0.625-in ODx 0.035-in wall tube. A man- 



(13) 


H D *! D o) 


These holes can be used to embed temperature probes in the wall of the inner test section tube. Five slots 


drel is slid inside a large bore sheath tube with a 0.75-in OD and 0.02-in wall. This assembly is drawn 
until the outer tube touches the inner tube. The assembly is then drawn again to the desired OD of the 
sheath tube. The outer surface can then be machined or otherwise finished. 
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Draw 



Figure 16. Inner tube assembly: (a) Predrawn and (b) postdrawn. 


Test section wall temperatures can be measured by inserting sensors into the axial holes of the 
inner tube. Two sensor types are possible. Ungrounded type K thermocouples with special error limits 
having 0.01-in OD sheaths, with a single bead each, allow for up to 10 wall temperature measurements 
to be made if thermocouples are fed from both ends of the inner tube. Fiber optic temperature sensors, 
0.005-in diameter, contained inside a 0.01 -in-diameter stainless steel capillary tube can also be used. 
These sensors, marketed by Luna Innovations, exploit Fabry-Perot interferometry to measure tempera- 
ture to within ±0.5 °C up to 1,200 °C. Temperature is measured at gratings spaced along the fiber at 
1-cm intervals, enabling a single fiber, fed from one end of the inner tube, to make measurements at up 
to 50 locations. If fiber optic sensors are placed in all five axial holes of the inner tube, up to 250 temper- 
ature measurements are possible. Use of thermocouples or fibers in an electrically heated tube requires 
isolation of the metallic sheath from the grounded jacket and manifold. 

Surface temperatures must be inferred from the temperature sensor embedded in the wall of the 
electrically heated inner tube. To get an accurate reading, the measurement must be corrected for the 
temperature drop between the thermocouple and the surface. Assume the region around the temperature 
sensor and the cooled surface can be approximated as a plane wall with uniform energy generation 
per unit volume, adiabatic on one face and convectively cooled on another. The surface temperature 
is related to the wall temperature by 


Tw=T, 
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(14) 


where the inner tube wall thickness is R and the distance of the thermocouple from the adiabatic face 
is r. Given the proposed test conditions, it is anticipated that this correction will be ~2 °C, for an overall 
temperature difference between the surface and the gas of 30 °C. 


For a constant surface heat flux condition, mean gas temperature as it flows through a channel 
is a linear function of position: 


Tg i x ) ~ + 


qnd 0 


men 


-x . 


(15) 
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If the outlet temperature and channel length are known, the relation between temperature and position 
can be simplified to 


T g (x) = T i + (T 0 -T i )j-, 0<x<L . (16) 

Provided power is uniformly distributed in the wall, this method of gas temperature calculation 
should prove satisfactory, with no disturbances to the flow field or ambiguity in the measurement of 
stream temperature. A contingency, should this approach not work, is to place thermocouples into the 
pressure taps so that they are flush mounted and out of thermal contact with the jacket. 
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5. DATA REDUCTION 


5.1 Pressure Drop 

The annuli tested in this study will have a finish better than 1 pm, yielding el d= lx 10 -4 , which 
is effectively smooth at Re = 30,000. For laminar flow, the entrance region is taken to be L e /d~ 0.06 Re. 
The maximum laminar entrance length, at a Reynolds number of 2,300, is 138 diameters, the longest 
length possible regardless of regime. For fully turbulent flow, the hydrodynamic entrance region is 
L e ld~4.4 Re 1/6 . Table 4 compares the entrance length to hydraulic diameter ratios for the four annular 
test sections in the 4,000 to 100,000 Reynolds number range using data recommended for annular flow 
by Jonsson and Sparrow. 23 


Table 4. Flydrodynamic entrance lengths for test sections 
in turbulent flow. 


Lld h 

Re 

144 

AF1 

95 

AF2 

72 

AF3 

29 

AF4 

4,000 

22 

19 

18 

15 

10,000 

25 

22 

21 

18 

30,000 

30 

27 

25 

21 

100,000 

37 

33 

31 

26 


Pressure measurements will be made at 2-in intervals along the 18-in length of the test section. 

In the laminar regime, flow will be developing through the entire length of the test section and a fully 
developed friction factor will not be reached. Above the transition point when Re> 2,300, with hydraulic 
diameter <0.25 in, fully developed flow exists along a significant fraction of the test section length. 

The test section with the largest hydraulic diameter, AF4, has fully developed flow over less than half 
its length. It is expected that pressure drop across the annular test section will be fully accountable 
by hydraulic radius correlation once gas viscosity and Reynolds number at composition is properly 
established. 


Kays uses a relation for the friction factor along the fully developed length of a heated annular 
channel. 24 This relation accounts for friction- and temperature -related acceleration effects: 
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(17) 


For flows with measurable heat transfer, the mass flow rate can be inferred from the electrical measure- 
ments and temperature rise across the test section or the volumetric flow meter. For adiabatic flows, 
mass flow rate can only be established with the flow meter. Estimates of uncertainty intervals for gas 
friction factor can be found in appendix C. 
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5.2 Test Section Thermal Equilibrium 

Thermal equilibrium in the test article occurs when thermal gradients across the test section no 
longer have a time-varying component. A one-dimensional transient conduction analysis for the jacket 
tube and end flange provides an order of magnitude estimate of the time required for the test section 
to reach equilibrium after a step change in boundary condition. The Biot number, Bi = hLk~ l , for the 
jacket tube was calculated. Here, the jacket wall thickness is used for the length scale, L, in the Biot 
number of the jacket. The radial distance between the ID of the jacket tube and the OD of the flange 
is the length scale of the end flanges. Heisler charts 25 were used to find the internal energy change as 
a function of time for a plane wall of thickness 2 L. This yields the Fourier number corresponding to 
an internal energy change of Q/ Q n = 0.99. The time constant for this change was then calculated using 
T = p s c s k s R Fo. Table 5 shows the results for the jacket including Biot numbers, Fourier numbers, 
conduction distances, and resulting time constants for the annulus ratios under consideration. 


Table 5. Jacket time constants. 


Variable 

AF1 

AF2 

AF3 

AF4 

Units 

Description 

Ps 

7,900 

7,900 

7,900 

7,900 

kg-m -3 

Density of stainless steel 

Cs 

400 

400 

400 

400 

J/kg-K 

Heat capacity of stainless steel 

k s 

17 

17 

17 

17 

W/m-K 

Thermal conductivity of stainless steel 

h 

1,245 

819 

623 

249 

W/m-K 

Gas heat transfer coefficient 

dr 

0.25 

0.25 

0.25 

0.25 

in 

Jacket wall thickness 

Bi 

0.465 

0.306 

0.233 

0.93 

- 

Biot number for jacket 

Fo 

18.5 

22.4 

26.4 

30.3 

- 

Fourier number for jacket 

T 

139 

168 

208 

779 

s 

Time constant for jacket 


Table 6 shows similar results for the flange. The time constant for the flange is about an order of 
magnitude higher than for the jacket. This reflects the larger volume of the flange relative to the surface 
area available for heat removal. The jacket tube time constant more accurately describes steady state for 
the test section. Use of the longer duration flange time constant may be considered a sufficient condition 
for thermal equilibrium. 


Table 6. Flange time constants. 


Variable 

AF1 

AF2 

AF3 

AF4 

Units 

Description 

Ps 

7,900 

7,900 

7,900 

7,900 

kg-m -3 

Density of stainless steel 

Cs 

400 

400 

400 

400 

J/kg-K 

Heat capacity of stainless steel 

k s 

17 

17 

17 

17 

W/m-K 

Thermal conductivity of stainless steel 

h 

1,245 

819 

623 

249 

W/m-K 

Gas heat transfer coefficient 

dR 

1.625 

1.593 

1.563 

1.375 

in 

Radial distance jacket OD to flange 

Bi 

3.024 

1.949 

1.454 

0.512 

- 

Biot number for flange 

Fo 

4.4 

7.6 

9 

15.3 

- 

Fourier number for flange 

T 

1,386 

2,321 

2,632 

3,464 

s 

Time constant for flange 


26 


5.3 Heat Transfer 


The test section is assumed sufficiently insulated to keep heat loss a small fraction of the applied 
power. An electrically heated tube provides power to the test section for the single-channel heat transfer 
test. This power is determined using measurement of voltage across the test section terminals and the 
voltage across a known resistance in the circuit. Using heating rate established in equation (13), the cor- 
rection for temperature drop across the tube wall in equation (14), and the computed gas temperature in 
equation (15), the Nusselt number is 


Nu = 


A ( T w ~ T g ) k g 


( 18 ) 


where k g is the thermal conductivity of the gas, d h is the hydraulic diameter of the test section channel, 
q/A is the power supplied per unit area supplied to the gas locally, T is the measured wall temperature, 
and T is the gas bulk temperature. Gas thermal conductivity is assessed based on methods outlined 

o 

earlier in this Technical Memorandum (TM) at the bulk temperature. Natural convection effects in the 
annular channel should be quite small, regardless of orientation, if the jacket tube is well insulated. Esti- 
mates of uncertainty intervals for Nussult number at various conditions can be found in appendix C. 
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6. HARDWARE PROCUREMENTS 


For all hardware acquisitions that exceeded a cost of $2,500, the standard full and open competi- 
tive procurement process was implemented. Once final vendor submissions were received, selection was 
based on evaluation criteria consisting of a combined cost and technical evaluation basis. Specifications 
for many of the hardware components have also been developed and procurements initiated. The larger 
components that were procured include the following: 

• Laminar flow meters. 

• Gas chiller. 

• Test section interior tube. 

• Test section jacket. 

• Inert gas system. 

• Gas compressor or booster pumps. 

• Direct current power supplies. 

Smaller purchases <$2,500 for items such as tubing, relief valves, fittings, raw materials, vacuum 
system components, etc., have also been made and are not listed individually. 
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7. STATUS AT CLOSEOUT AND CONCLUSIONS 


This TM has provided an overview and critical evaluation of the test design to date. Conven- 
tional relations found in the literature predict that fully developed conditions for turbulent flow should 
be established over a sizable fraction of the test sections with annulus ratios >0.7. Entrance region 
effects for turbulent flow will be an important consideration in the interpretation of data for a test section 
with an annulus ratio of 0.5. Time constants to test section equilibrium should normally be on the order 
of an hour after a large change in boundary conditions. Several important findings are included in the 
error analysis contained in appendix C. Below annulus ratios of ~0.8 error for measured Nusselt number 
should be less than ±9%. Due to constraints on fabrication tolerances, test section with annulus ratios 
>0.8 are unlikely to yield well-controlled measurements of Nusselt number. Uncertainty intervals for 
Nusselt number are expected to correlate only weakly with Reynolds number. 

This work might be made to cover a wider range of test conditions and geometrical configura- 
tions. These extensions might include different gas mixtures, temperatures, pressures, mass flow rates, 
annulus ratios, and heat transfer enhancement methods. Every attempt was made in the initial design to 
permit later accommodation of these contingencies with minimal changes to the basic hardware configu- 
ration. In its current embodiment, the apparatus can make accurate heat transfer measurements for circu- 
lar annuli with an internally heated tube having 0.1 <e<0.8. Fabrication tolerances tighter than ASTM 
standards are required to obtain good data for annuli with e>0.8. Tests of circular tubes (e=0) or circu- 
lar annuli that are heated on the inner and outer surfaces are also possible by direct heating of the jacket 
tube. With the existing glass-to-metal power feeds, the test section should be operable near room tem- 
perature to at least 250 psia. By active cooling of the feeds and using all welded construction, the current 
stainless steel test section could be brought to at least 700 °C. Operation at these temperatures requires 
modification to the inner tube power lead located at the exit of the annulus. For testing He-Ar mixtures, 
where some gas wastage is tolerable, the RIX compressor should be satisfactory. Tests with He-Xe mix- 
tures will require a booster pump or hermetically sealed turbo blower with very low leakage rates. Cost 
constraints could justify the building of a liquid nitrogen-cooled trap to recover Xe from the flow loop. 

Although this project was terminated prior to setting up the hardware support systems within the 
laboratory, two of the basic support systems were in fabrication: (1) The test section and (2) the gas cir- 
culation loop, both of which may be used to support the future evaluation of heat transfer in gas-cooled 
core channels. Four test sections for single-channel testing have been designed and one was built using 
MSFC-directed funds. Engineering of the leak-tight compressor is still required before testing with He- 
Xe in the loop. Also, design of the connection of the power supplies to the test section still needs to be 
finalized. The power supplies are currently well within their rated limits. Outside of the length of the gas 
circulation loop lines, there should be few constraints on the layout of the system. 
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APPENDIX A -THERMOPHYSICAL PROPERTY FITS 


Third-order polynomial fits for thermal conductivity and viscosity of Ar gas with deviations from 
accepted values over the range 300 K < T < 2,000 K are as follows: 


/c Ar (r) = 0.0011434 + 6.1189xl0“ 5 r-2.5315xl0“ 8 r 2 + 5.9271xl0“ 12 r 3 , (19) 

with error limits 1.959/-0.414%, 
and 

jl Ar (T) = 1.5353 x 10 -6 + 7.8135 x 10“ 8 T - 2.2451 x 10 _11 r 2 + 6.643 x 10“ 15 T 3 , (20) 

with error limits 2.137 1-0.611%. 

Figures 17 and 18 compare curve fits of the thermal conductivity and dynamic viscosity of Ar 
as a function of temperature between 100 and 2,000 K to accepted values found in Touloukian. In both 
cases, third-order polynomial fits were sufficient to capture the data with a maximum deviation of 2%. 



Figure 17. Thermal conductivity of Ar gas 
at low density: accepted values 
from Touloukian with third-order 
polynomial fit. 3 



Figure 18. Dynamic viscosity of Ar gas 

at low density: accepted values 
from Touloukian with third-order 
polynomial fit. 4 
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Third-order polynomial fits for thermal conductivity and viscosity of He gas with deviations 
from accepted values over the range 300 K < T < 2,000 K with units W/m-K and kg/m-s, respectively, 
are as follows: 

k He (r) = 0.036953 + 3.9178 xl0“ 4 r- 8.1751 xl0“ 8 r 2 + 1.5913x 10“ n r 3 , (; 

with error limits 1 .562 /—1 .416%, 
and 

^ He (r) = 4.693 xl0“ 6 + 5.463 xl0“ 8 r- 1.7247 xlO _11 T 2 + 3.159 xl0“ 15 T 3 , (; 

with error limits 1.382/-0.434%. 

These fits are compared with accepted values in figures 19 and 20. 
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Figure 19. Thermal conductivity of He gas 
at low density: accepted values 
from Touloukian with third-order 
polynomial fit. 3 


Figure 20. Dynamic viscosity of He gas 

at low density: accepted values 
from Touloukian with third-order 
polynomial fit. 4 


Third-order polynomial fits for thermal conductivity and viscosity of N 9 gas with deviations 
from accepted values over the range 300 K< T< 2,000 K are 

% 2 (T) = 2.7679 x 10“ 4 +9.8753x10“ 5 T- 5.1554 x10“ 8 T 2 +1.5427 xlO _11 r 3 , (23) 

with error limits 1.157 / -0.434%, 
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and 


ju Ni (t) = 1.9334 x 10“ 6 + 5.8637 x 10“ 8 T - 2.5153 x 10 _11 r 2 + 5.2386 x 10“ 15 T 3 , (24) 

with error limits 2.563 /-0. 771%. 

Again, the units on temperature are K, and conductivity and viscosity are W/m-K and kg/m-s, 
respectively. Figures 21 and 22 plot these fits compared to accepted values as a function of temperature. 
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Figure 21. Thermal conductivity of N 0 gas 
at low density: accepted values 
from Touloukian with third-order 
polynomial fit. 3 


Figure 22. Dynamic viscosity of N 0 gas 
at low density: accepted values 
from Touloukian with third-order 
polynomial fit. 4 


Third-order polynomial fits for thermal conductivity and viscosity of Xe gas with deviations from 
accepted values over the range 300 K < T < 2,000 K: 

k Xe {T ) = 1.1526 x 10“ 4 + 1.9695 x 10“ 5 T - 4.8985 x 10“ 9 r 2 + 6.5918 x 10“ 13 T 3 , (25) 

with error limits 1.157 /-0.434%, 
and 

/l Xe ( T ) = -9.921 x 10~ 7 + 9.0488 x 10“ 8 T - 3.3146 x 10 _11 r 2 + 6.8634 x 10“ 15 T 3 , 
with error limits 2.563/-0.771%. 

These fits are compared to accepted values found in the literature in figures 23 and 24. 
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Figure 23. Thermal conductivity of Xe gas 
at low density: accepted values 
from Touloukian, Jain, and Jody 
with third-order polynomial fit. 3 ’ 26 ’ 27 


Figure 24. Dynamic viscosity of Xe gas 

at low density: accepted values 
from Touloukian with third-order 
polynomial fit. 4 


Figures 25 and 26 plot deviations for each polynomial curve fit from accepted value as a function 
of temperature. Above 300 K, deviations in viscosity are generally <1%. Deviations in conductivity are 
somewhat large, yet above room temperature, stay within a 1.5% band. These deviations are within the 
test uncertainty cited for the thermal conductivity of gas found in a survey by Ho et al. 28 



Figure 25. Deviation in third-order thermal 
conductivity curve fits from 
accepted values found in 
Touloukian and Jody. 3 ’ 27 



Figure 26. Deviation in third-order dynamic 
viscosity curve fits from accepted 
values found in Touloukian 
and Jody. 4 ’ 27 
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APPENDIX B- DRAWINGS OF TEST SECTION COMPONENTS 


Figures 27-40 show various schematics of the test section components. 
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Figure 27. Single-channel test section assembly. 
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Figure 28. Inlet flange manifold. 
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Figure 29. Inlet flange manifold— pipe with spider attachment. 
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Figure 31. Spider attachment. 
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Figure 33. Inlet flange. 
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Figure 34. Inner tube plug. 
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Figure 35. Jacket pipe. 
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Figure 36. Inner tube predraw assembly. 


44 


1 


u 


CnI 


-fr 





Figure 37. Inner tube postdraw assembly. 
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Figure 38. Exit manifold. 
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Figure 39. Outlet centering ring. 


47 



Figure 40. Outlet manifold. 
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APPENDIX C-TEST UNCERTAINTY ANALYSIS 


Given the restrictions on resources, this test series will be predominately single sample with few 
replicates. In this section, the uncertainty intervals for key results are estimated using the method to 
describe propagation of uncertainties in single-sample experiments found in Kline and McClintock. 29 
Uncertainty intervals are estimated for the independent dimensionless groups Reynolds and Prandtl 
number and the dependent dimensionless groups they correlate, friction factor, and Nusselt number. 
Here, all variables that comprise these groups are assumed normally distributed with 20 to 1 odds that 
the measurement falls within the interval about the mean. Steady-state conditions are assumed for the 
fully developed region near midpoint of the test section. Uncertainty intervals are calculated for each 
test section operating near room temperature. The sensitivity of these estimated uncertainty intervals to 
variations in test conditions, such as temperature and Reynolds number, is also examined. 

C.l Prandtl Number Uncertainty 

The Prandtl number, Pr=/ic p /k, of the gas is used here as an independent variable in the correla- 
tion of the nondimensional heat transfer rate. It is a function of specific heat— a thermodynamic property, 
and also viscosity and thermal conductivity— both transport properties. The shape of the Prandtl num- 
ber versus composition curve for He-Ar mixtures is shown in the introductory paragraphs of this TM. 
Prandtl number is generally a weak function of temperature and, near its minima, a weak function of 
composition. Given these relations, the variation of Prandtl number should be insensitive to small varia- 
tions in temperature and composition. The same does not hold for the variables that comprise Prandtl 
number. Gas viscosity and thermal conductivity depend on temperature and composition. Ideal gas 
specific heat depends on composition only. Viscosity and thermal conductivity correlations agree with 
data for He-Ar gas mixtures reported to the 1% and 4% levels, respectively, over the range of conditions 
considered here. Property correlations using such published values should be accurate to these levels for 
a given temperature or composition. The associated uncertainty interval is 


8Pr 

Pr 


\2 ( 


8/i 

P 


V ( 
+ 




8c, 


V 


c„ 

\ P J 


+ = 0.01 2 +0.01 2 +0.04^ 


= 0.0424 , 


(27) 


where the uncertainty intervals for viscosity and temperature are estimated based on agreement between 
correlation and data. The uncertainty interval for specific heat is estimated based on possible gas compo- 
sition errors. To reduce the chance of inadvertent contamination, the test section will be evacuated 
<10 -6 torr before premixed, high-pressure gas is added. The composition of bottles containing He-Ar 
mixtures is typically known to the part-per-million level. For this estimate, gas composition will be 
assumed known to within 1%. This rough estimate yields a value of Prandtl number near the He-Ar 
compositional minima of Pr=0 . 39 + 0.017 (20 to 1 odds). 
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C.2 Friction Factor Uncertainty for Adiabatic Flows 

The friction factor for fully developed turbulent flow through an annulus is a weak function of 
eccentricity when c^-^O. 8 Since it is possible that friction factor measurements will not be made on a 
heated test section, applied power and temperature rise through the test section cannot be relied upon to 
establish the flow condition. Instead, a laminar flow meter provides the mass flow rate measurement for 
the adiabatic case. The friction factor for adiabatic flows between two points along the length of the test 
section can be described as f- 2Ap(Ax)~ 1 A 2 V~ 2 p ~ 1 d h . The resulting uncertainty interval for friction factor 
in terms of these variables is: 



Pressure taps are placed at eight locations along the test section at 5.08-cm intervals. A Pressure 
Systems Esterline model 9116 network analyzer senses differential pressure at each tap with respect to 
a reference port located near the center of the test section. The manufacturer quotes a measurement reso- 
lution of +0.003% full scale, a static accuracy of +0.15% full scale for p < 20 kPa, and ±0.05% full scale 
for p > 20 kPa. Total thermal error is ±0.0015% full scale per degrees Celsius for p < 20 kPa and ±0.001% 
full scale per degrees Celsius for p > 20 kPa. The line pressure effect for this sensor is ±0.0003% full 
scale per psi for absolute p> 100 kPa. If the pressure drop between taps is assumed <20 kPa and thermal 
error is taken to be 20 °C, the relative uncertainty for a 7-kPa, full-scale sensor measuring a 224 Pa 
reading is 


8Ap 7,000 Pa 
A p 224 Pa 


22»V + f^ET 20 °c 

100 ) l 100 °c 


+ 


0.0003 

100 


— lOOpsi I 
psi ) 


-i 0.5 


0.0496 . (29) 


The uncertainty interval for the 2-in spacing between pressure taps is taken from tolerances on the 
fabrication drawings: 


8Ax 

Ax 


0.010 

2 


0.005 . 


(30) 


The uncertainty interval for the channel cross-sectional area, A = k[d} - <7 2 ] / 4 are taken from 
the fabrication drawings: 


SA 

A 


'8DA 

{ D i J 




0.002^j 2 C 0.005 ^ 2 
0.815 J + 1 0-625 J 


0.00886 . 


(31) 


Table 7 shows estimated uncertainty for adiabatic friction factors for He-40%Ar near room tem- 
perature based on an assumed Blasius-type friction factor correlation. The second column in the table 
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contains the values used to compute uncertainty in the column labeled AF2. The estimates for AF1, 

AF2, and AF3 use Esterline sensors with 7-kPa full scale. For AF4, a 2-kPa, full-scale sensor is used. 
For low annulus ratio, AF4, the pressure drop is quite low and the sensitivity of the pressure transducer 
dominates uncertainties that approach 30%. As the annulus ratio increases, uncertainties in gas density 
and hydraulic diameter become important. The uncertainty for these remaining cases, with 20 to 1 odds, 
is under 6%. Table 8 shows the influence of Reynolds number on uncertainty in pressure drop measure- 
ment. More sensitive pressure transducers will be needed if accurate readings are desired at Re< 10,000. 


Table 7. Uncertainties in friction factor measurement for r* =0.76, He-40%Ar. 


Symbol 

Value 

Unit 

Description 

AF2 

P 

689,303 

Pa 

Absolute pressure 

0.0007 

Ap 

303 

Pa 

Pressure drop 

0.0367 

Ax 

2 

in 

Spacing between pressure taps 

0.005 

Ps 

5.237 

kg/m 3 

Gas density 

0.0203 

m 

0.0144 

kg/s 

Mass flow rate 

0.0086 

A 

0.2148 

in 2 

Cross-sectional area 

0.0084 

d h 

0.19 

in 

Hydraulic diameter 

0.0186 

Re 

30,000 

- 

Reynolds number 

0.0431 

f 

0.007 

- 

Friction factor correlation 

0.052 


Table 8. Uncertainties in friction factor measurement for He-40%Ar mixtures. 


Re (?) 

0.83 

0.76 

0.71 

0.5 

3,000 

0.344 

1.201 

2.73 

15.2 

7,000 

0.089 

0.275 

0.622 

3.46 

15,000 

0.047 

0.089 

0.167 

0.91 

22,000 

0.044 

0.052 

0.091 

0.47 

30,000 

0.043 

0.052 

0.06 

0.27 


C.3 Nusselt Number Uncertainty 

The test section is assumed sufficiently insulated to keep heat loss a small fraction of the applied 
power. An electrically heated tube provides power to the test section for the single-channel heat transfer 
test. This power is determined using measurement of voltage across the test section terminals and the 
voltage across a known resistance in the circuit. From Rohsenow, Hartnett, and Ganic, 8 a simple linear 
relation between Nusselt number and eccentricity can be assumed, provided the annulus ratio is <0.8. 
Expressing Nusselt number in this way produces the following formulation: 


qd h (\-e/2) 
A (t w ~ T g ) k g 


(32) 


where k g is the thermal conductivity of the gas, d h is the hydraulic diameter of the test section, q/A is the 
power supplied per unit area supplied to the gas locally, T w is the measured wall temperature, T is the 
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gas bulk temperature, and e is the eccentricity of the annulus. The temperature difference is based on the 
reading of wall thermocouples, with accuracy on the order of ±1.1 °C. Gas properties are assessed based 
on methods outlined earlier in this TM at the local temperature. Natural convection effects in the annular 
channel should be quite small, regardless of orientation, since the jacket tube is well insulated. 

C.3.1 General Uncertainty Equations 

For the current uncertainty analysis, the expression in equation (32) is used and the Nusselt 
number can be expressed as a function of the following terms: 


Nu = f(q,A,d h ,T w ,T g ,k g ,e) . 
The linear propagation of errors results in 


(33) 


CA7 dNu c dNu c dNu c , dNu Qrn dNu dNu dNu c 

SNu = 8q + 8 A + 8cl h + 8T W + 8T„ + 8k„ + 8e 

dq 8A dd h h dT w w dT„ 8 })k„ 8 de 


(34) 


Squaring, expanding, and setting all correlated terms to zero results in the final expression: 




1 dNu 2 

v^J 5r ”' 


+ 


1 dNu} 2 


ST} + 


K dk g ) 8 \ 8e ) 


(35) 


All partial derivatives listed in equation (35) are as follows: 


dNu _ 


e/2) 

dq 

a( T w - 

T g) k g 

dNu _ 

qd h ( 1 

- e/2 ) 

dA 

a 2 (t„- 

~ T g) k g 

dNu 

q{ i- 

e/2) 

ddh 

a(t w - 

T g) k g 

dNu _ 

qd h (l 

-e/2) 


dT, 


(36) 


(37) 


(38) 


2 


(39) 


( 40 ) 


and 


dNu 

qd h (\-e/2) 

dT g ~ 

A ( T w~ T g) k g 

dNu _ 

qd h (l-e/2) 

dk g 

A ( T « ~ T s) k l 

dNu 

1 dd/i 

de ' 

2A(T W -T g )k g 


(41) 


(42) 


Substituting equations (36)-(42) into equation (35) and rearranging each derivative term by 
dividing it by its square results in the following: 


8Nu 2 


qd h (l - e/2) 

< A { T w~ T g) k g , 




N 9 ( \ 

SA\ 2 + - 8T W 

A> 


^ ) 

2 

+ 


(t -T ) 

k„ 

\\ w s)) 


\ 8 J 


-8e 


2(1- e/2) 


(43) 


Each of the terms in equation (43) must now be determined. 

C.3.2 Heating Rate Uncertainty Interval 


The heating rate is established assuming electrical power is uniformly dissipated to the inner tube 
with homogeneous composition and cross section. Additional heating or cooling can occur from para- 
sitic losses across insulation that couples the test section gas and the environment: 


2 jckgLlTg-Tj) 


In (D s /D a ) 


(44) 


A layer of 0.068-W/m-K insulation ~2.3 in thick on a 1 -in-diameter tube is sufficient to keep 
the second parasitic term to <4 W. The uncertainty interval is 


8q 2 = 


N) se 2 + 

dE ) 


g K * 8 - 1 


(45) 
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The nominal current through the test section is 1,500 A + (0. 13/100 1,500) A as measured with a high- 
precision resistor. The terminal voltage is 0.5 V ± (0.03/100 0.5) V. The thermal loss from the test 
section assumes use of vermiculite, or similar packed or loose-packed insulation, and an overall temper- 
ature gradient between the gas and ambient established at the entrance of the test section: 


(46) 

(47) 

C.3.3 Hydraulic Diameter Uncertainty Interval 

The hydraulic diameter, d, uncertainty interval is estimated using fabrication tolerance: 


2 

8q 2 = I 2 8E 2 + E 2 8l 2 + 8q 2 = (l, 500 x 0.15 x 10 -4 j~ +(0.5 x 1.95) 2 + 3.7 2 
or 

g = 888 + 4.05 W (20 to 1) . 


dfo — Di — d a , (48) 

(8d h ) 2 =(8Di) 2 +(8d 0 ) 2 =(0.005) 2 +(0.002) 2 , (49 ) 

or 

d = 0.19 ±0.00539 in (20 to 1) . (50) 

C.3.4 Heat Transfer Area Uncertainty Interval 

Power is generated at a uniform rate in the inner tube by electrical dissipation. The convective 
heat flux is based on the wetted area of the inner tube given by A=nd a L, where the length is along 
the test section enclosed by the jacket: 


or 



(51) 


A = 35.3 ±0.344 in 2 (20 to 1) . 


(52) 


C.3.5 Surface Temperature Uncertainty Interval 

Surface temperatures are measured with a thermocouple embedded into the wall of the electri- 
cally heated inner tube. To get an accurate reading, the measurement must be corrected for the tempera- 
ture drop between the thermocouple and the surface, assuming the region around the thermocouple and 
the cooled surface can be approximated as a plane wall with uniform energy generation per unit volume, 
adiabatic on one face and convectively cooled on another. The surface temperature is related to the 
thermocouple temperature by 
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(53) 


T w =T t 


mrt 

q R 


2 ( 


2 K 


R 


0 <r<R 


where the inner tube wall thickness is R and the distance of the thermocouple from the adiabatic face 
is r. Given the test conditions, it is anticipated that this correction will be ~2 °C for an overall tempera- 
ture difference between the surface and the gas of 30 °C. The uncertainty interval for the wall thermo- 
couple is 


(«’*) 2 = + 


dr. 


dT u 


dT u 


dT u 


n 8q"'\ + + ^8k s I +1 n 

dq'" 1 ) { dR ) [ dk, s ) [dr 


8r 


The partial derivatives for each term are as follows: 
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dT t 


JT W 
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v 0.084 j 
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~YT ////}2 ( 2 ^ 

dT w _ q R , r 


dk s 2 k 2 


R z 


(4)888(0.084(0.0254)) 
tt(o. 625 2 - 0.457 2 )l8(0.0254) 3 17 

(4)888 (0.084(0.0254)) 2 


7T ( 0.625 2 - 0.457 2 18(0.0254) 3 2(17) 2 


^ 0.049 2 ^ 

v 0.084 2 j 


dT w _ q"'r _ (4)888(0.049(0.0254)) 

-fr-~k7- Ja,oc2 .k, 2\ 10 ,nnoc^3 1 


and 

dT w = q"'r 

k s 7T (^0.625 * - 0.457 ^ 118(0.0254) J ll 
The uncertainty interval for the thermocouple measurement is 

(8T t ) 2 = (8T Wire ) 2 + (<5r IcePoint ) 2 + (<)T RoLincl ) “ + (dT Attach ) 2 , 
(« 8T t ) 2 = (1.10) 2 + (0.05) 2 + (0.01) 2 + (0.10) 2 , 


(54) 


(55) 


(56) 


(57) 


(58) 


(59) 


(60) 

( 61 ) 
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and 


r,=321.4±l.ll K(20 to 1). 


(62) 


The uncertainty intervals for the tube wall thickness, R, and radial position or the thermocouple, r, 
are taken from fabrication drawings: 


{SR) 2 = (0.015 / 2) 2 + (0.005 / 2) 2 and (Sr) 2 = (0.015 / 2) 2 . (63) 

The uncertainty interval for thermal conductivity of the solid tube wall is assumed negligible. 

C.3.6 Gas Temperature Uncertainty Interval 


For a constant surface heat flux condition, mean gas temperature as it flows through a channel 

/ \ qnd n 

is a linear function of position, T g (x) = 7) + — x. If the outlet temperature and channel length 

are known, the relation between temperature and position can be simplified to 

T g (x) = T i +(T„-T i )j , (64) 

where 0 <x< L. 

The uncertainty interval is then 



Substituting values into this equation yields 



7^ = 291.1 + 1.12 K (20 to 1). (67) 


C.3.7 Gas Thermal Conductivity Uncertainty Interval 

Values for the thermal conductivity of gas mixtures have not been firmly established over a wide 
temperature range by either theory or test. At low pressure, gas thermal conductivity is a function of 
its composition and temperature. Deviations between measured and theoretical values can often range 
+10%. For this estimate, the functional relation described in appendix A will be assumed to correspond 
to the data and an uncertainty interval will be established based on uncertainties in temperature and com- 
position using numerically calculated partial derivatives: 
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( 68 ) 
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A 8T, 


dr, 
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d _K 

dn 


-8n 


2 _[ (0.0648-0.0584) 
dkg ~{ (311-271) 


2 

+ 


(0.0754-0.0505) 
(0.3-0. 5) 



2 


and 

k =0.06170±0.00249 W/m-K (20 to 1) . 

o 

C.3.8 Annulus Eccentricity Uncertainty 


(69) 


(70) 


Annulus eccentricity is the distance between the centerlines of the inner and jacket tubes divided 
by the mean gap between the tubes. Eccentricity can affect local Nusselt numbers as shown in Rohse- 
now, Hartnett, and Ganic. 8 A linear relation is adopted between Nusselt number and eccentricity. Eccen- 
tricity is estimated based on measurements of total indicated runout of the inner and jacket tubes which 
are 0.0085 in and 0.005 in, respectively: 


2e 

D ; - dr 


2 [{ tir i ) 2 +( tir o) 

D: - dr. 


1/2 


2 ((0.0085) 2 +(0.005) : 
0.815-0.625 


1/2 


(71) 


e = 0 + 0.104 (20 to 1) . 


(72) 


C.3.9 Reynolds Number Uncertainty 

The Reynolds number for a fluid flowing through a circular annulus, Re = 4 q [ cj T () - 7)) n{D t 
+ d 0 )n] L , can be described as a function of heating rate, temperature rise across the test section, the 
principle diameters of the passage, gas specific heat, and viscosity. It is used here as an independent vari- 
able in correlation of the nondimensional friction factor and the heat transfer rate. 

Uncertainty in experimental variables and results at r* = 0.833, r*=0.878, r* = 0.714, and 
r* = 0.500 are shown in tables 9-12, respectively. 
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Table 9. Uncertainty in experimental variables and results at r* = 0.833. 


Variable 

Value ± Uncertainty 

Description 

Interval 

Gas 

He-40Ar ± 2% 

Gas composition 

0.02 

E 

2.0518 ± 6.16X10 -4 V 

Terminal voltage 

0.0003 

1 

413 ±0.54 A 

Current through test section 

0.0013 

R 

4.97X10- 3 ± 6.63 x10“ 6 ohm 

Resistance of test section 

0.0013 

q 

848 ± 3.78 W 

Power supplied to inner tube 

0.0045 

q 0 

-3.61 ± 0.63 W 

Heat loss (insulation only) 

0.1751 

q" 

2.01 x10 7 ± 6.86 xIO 5 W/m 3 

Volumetric heating rate 

0.0341 

T t 

330.27 ±1.1 K 

Embedded thermocouple temperature 

0.0293 

T w 

328.49 ± 1.61 K 

Surface temperature at tube OD 

0.043 

T 0 

311 ± 1.1 K 

Outlet temperature from test section 

0.0293 

\ 

300 ± 1.1 K 

Room temperature 

0.0293 

T 9 

291 ±1.12 K 

Bulk gas temperature 

0.0299 

T, 

271 ± 1.1 K 

Inlet temperature to test section 

0.0293 

T W ~Tg 

37.49 ± 1.96 K 

Wall to gas temperature difference 

0.0523 

TfT w 

1.78 ± 1.18 K 

Thermocouple to surface drop 

0.6615 

k 

25 ±0.01 K 

Test section total length 

0.0004 

k 

18 ±0.01 K 

Test section active length 

0.0006 

X 

9 ± 0.1 in 

Measurement position 

0.0111 

Ds 

6 ± 0.1 in 

Insulation outside diameter 

0.0167 

D 0 

1.25 ±0.005 in 

Jacket outside diameter 

0.004 

Di 

0.75 ± 0.002 in 

Jacket inside diameter 

0.0027 

Di 

0.75 ± 0.005 in 

Jacket total indicated runout 

0.0067 

d o 

0.625 ± 0.005 in 

Inner tube outside diameter 

0.008 

d o 

0.625 ± 0.009 in 

Inner tube total indicated runout 

0.0136 

d 

0.555 ± 0.005 in 

Thermocouple diameteral location 

0.009 

d i 

0.457 ±0.015 in 

Inner tube inside diameter 

0.0328 

d h 

0.125 ±0.004 in 

Hydraulic diameter 

0.0282 

R 

0.084 ±0.016 in 

Inner tube wall thickness 

0.1882 

r 

0.049 ±0.016 in 

Thermocouple to adiabatic face distance 

0.3227 

A 

35.33 ± 0.344 in 2 

Heat transfer area 

0.0097 

k g 

0.066 ± 2.5 x 10- 3 W/m-K 

Gas thermal conductivity 

0.0379 

<9 

1,131 ±1.13x10 1 J/kg-K 

Gas specific heat 

0.01 

Pg 

2.28x10- 5 ±2.85x10“ 7 kg/m-s 

Gas viscosity 

0.0125 

c s 

400 ±4x10 J/kg-K 

Stainless steel heat capacity 

0.01 

k s 

17 ± 1.7x10 _1 W/m-K 

Stainless steel thermal conductivity 

0.01 

Ps 

7,900 ±7.9x10 1 kg/m 3 

Stainless steel density 

0.01 

r s 

7.2x10 -7 ± 7.2x10 -9 ohm-m 

Stainless steel resistivity 

0.01 

k i 

0.068 ±6.8x10“ 4 W/m-K 

Insulation conductivity vermiculite 

0.01 

e 

0 ± 0.158 

Eccentricity 

0.0789 

Re 

30,000 ±1,294 

Reynolds number 

0.0431 

Pr 

0.39 ± 0.02 

Prandtl number 

0.0404 

Nu 

47.8 ±91.1 

Nusselt number 

1.9059 
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Table 10. Uncertainty in experimental variables and results at r* =0.767 . 


Variable 

Value ± Uncertainty 

Description 

Interval 

Gas 

He-40Ar ± 2% 

Gas composition 

0.02 

E 

1.7817 ± 5.35x10^ V 

Terminal voltage 

0 

1 

498 ±0.65 A 

Current through test section 

0.001 

R 

3.58 x10“ 3 ± 4.77x1 0“ 6 ohm 

Resistance of test section 

0.001 

q 

888 ±3.91 W 

Power supplied to inner tube 

0.004 

% 

-3.73 ± 0.65 W 

Heat loss (insulation only) 

0.175 

q" 

2.11 xIO 7 ± 7.19x10 5 W/m 3 

Volumetric heating rate 

0.034 

T t 

352.54 ±1.1 K 

Embedded thermocouple temperature 

0.018 

T w 

350.68 ± 1.65 K 

Surface temperature at tube OD 

0.028 

T 0 

311 ±1.1 K 

Outlet temperature from test section 

0.018 

T a 

300 ±1.1 K 

Room temperature 

0.018 

T g 

291 ±1.12 K 

Bulk gas temperature 

0.019 

Ti 

271 ±1.1 K 

Inlet temperature to test section 

0.018 

Tw'Tg 

59.68 ± 2 K 

Wall to gas temperature difference 

0.033 

TfT w 

1.86 ± 1.23 K 

Thermocouple to surface drop 

0.662 

L, 

25 ±0.01 K 

Test section total length 

0 

La 

18 ±0.01 K 

Test section active length 

0.001 

X 

9 ±0.1 in 

Measurement position 

0.011 

D s 

6 ±0.1 in 

Insulation outside diameter 

0.017 

Do 

1.315 ±0.005 in 

Jacket outside diameter 

0.004 

Di 

0.81 5 ±0.002 in 

Jacket inside diameter 

0.002 

D, 

0.81 5 ±0.005 in 

Jacket total indicated runout 

0.006 

d o 

0.625 ± 0.005 in 

Inner tube outside diameter 

0.008 

d o 

0.625 ± 0.009 in 

Inner tube total indicated runout 

0.014 

d 

0.555 ± 0.005 in 

Thermocouple diameteral location 

0.009 

d i 

0.457 ±0.015 in 

Inner tube inside diameter 

0.033 

d h 

0.19 ± 0.004 in 

Hydraulic diameter 

0.019 

R 

0.084 ±0.016 in 

Inner tube wall thickness 

0.188 

r 

0.049 ±0.016 in 

Thermocouple to adiabatic face distance 

0.323 

A 

35.33 ± 0.344 in 2 

Heat transfer area 

0.010 

k g 

0.066 ± 2.5x10- 3 W/m-K 

Gas thermal conductivity 

0.038 

* 

1,131 ±1.13x10 1 J/kg-K 

Gas specific heat 

0.01 

Pg 

2.28x10- 5 ± 2.85 xIO” 7 kg/m-s 

Gas viscosity 

0.013 

C S 

400 ±4x10 J/kg-K 

Stainless steel heat capacity 

0.01 

k s 

17 ± 1.7x10 _1 W/m-K 

Stainless steel thermal conductivity 

0.01 

Ps 

7,900 ± 7.9 xIO 1 kg/m 3 

Stainless steel density 

0.01 

r s 

7.2 x 1 0 -7 ± 7.2 x 1 0 -9 ohm-m 

Stainless steel resistivity 

0.01 

k i 

0.068 ±6.8x10- 4 W/m-K 

Insulation conductivity vermiculite 

0.01 

e 

0 ± 0.104 

Eccentricity 

0.052 

Re 

30,000 ±1,293 

Reynolds number 

0.043 

Pr 

0.39 ±0.0166 

Prandtl number 

0.042 

Nu 

47.8 ±4.11 

Nusselt number 

0.086 
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Table 11. Uncertainty in experimental variables and results at r* =0.714. 


Variable 

Value ± Uncertainty 

Description 

Interval 

Gas 

He-40Ar ± 2% 

Gas composition 

0.02 

E 

1.8184 ± 5.46X10- 4 V 

Terminal voltage 

0 

1 

509 ±0.66 A 

Current through test section 

0.001 

R 

3.58 x 10“ 3 ± 4.77 x 10 -6 ohm 

Resistance of test section 

0.001 

q 

925 ± 4.04 W 

Power supplied to inner tube 

0.004 

q 0 

-3.84 ± 0.67 W 

Heat loss (insulation only) 

0.175 

q" 

2.2 x 10 7 ± 7.49 x10 5 W/m 3 

Volumetric heating rate 

0.034 

T t 

374.74 ±1.1 K 

Embedded thermocouple temperature 

0.013 

T w 

372.8 ± 1.69 K 

Surface temperature at tube OD 

0.021 

T 0 

311 ±1.1 K 

Outlet temperature from test section 

0.013 

T a 

300 ± 1.1 K 

Room temperature 

0.013 

T g 

291 ±1.12 K 

Bulk gas temperature 

0.014 

T, 

271 ± 1.1 K 

Inlet temperature to test section 

0.013 

Tw'Tg 

81.8 ± 2.03 K 

Wall to gas temperature difference 

0.025 

TfT w 

1.94 ± 1.28 K 

Thermocouple to surface drop 

0.662 

k 

25 ±0.01 K 

Test section total length 

0 

k 

18 ±0.01 K 

Test section active length 

0.001 

X 

9 ±0.1 in 

Measurement position 

0.011 

D s 

6 ±0.1 in 

Insulation outside diameter 

0.017 

D 0 

1.375 ±0.005 in 

Jacket outside diameter 

0.004 

Di 

0.875 ± 0.002 in 

Jacket inside diameter 

0.002 

Di 

0.875 ± 0.005 in 

Jacket total indicated runout 

0.006 

d o 

0.625 ± 0.005 in 

Inner tube outside diameter 

0.008 

d o 

0.625 ± 0.009 in 

Inner tube total indicated runout 

0.014 

d 

0.555 ± 0.005 in 

Thermocouple diameteral location 

0.009 

d i 

0.457 ±0.015 in 

Inner tube inside diameter 

0.033 

d h 

0.25 ± 0.004 in 

Hydraulic diameter 

0.014 

R 

0.084 ±0.016 in 

Inner tube wall thickness 

0.188 

r 

0.049 ±0.016 in 

Thermocouple to adiabatic face distance 

0.323 

A 

35.33 ± 0.344 in 2 

Heat transfer area 

0.01 

k a 

0.066 ±2.5x10- 3 W/m-K 

Gas thermal conductivity 

0.038 

C 9 

1,131 ±1.13x10 1 J/kg-K 

Gas specific heat 

0.01 

Pg 

2.28x10- 5 ±2.85x10- 7 kg/m-s 

Gas viscosity 

0.013 

C S 

400 ±4x10 J/kg-K 

Stainless steel heat capacity 

0.01 

k s 

17 ± 1.7x1 0 _1 W/m-K 

Stainless steel thermal conductivity 

0.01 

Ps 

7,900 ±7.9x10 1 kg/m 3 

Stainless steel density 

0.01 

r s 

7.2x10“ 7 ± 7.2x10- 9 ohm-m 

Stainless steel resistivity 

0.01 

k i 

0.068 ± 6.8 x10“ 4 W/m-K 

Insulation conductivity vermiculite 

0.01 

e 

0 ± 0.079 

Eccentricity 

0.039 

Re 

30,000 ±1,293 

Reynolds number 

0.043 

Pr 

0.39 ± 0.02 

Prandtl number 

0.04 

Nu 

47.8 ± 2.99 

Nusselt number 

0.063 
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Table 12. Uncertainty in experimental variables and results at r* =0.500. 


Variable 

Value ± Uncertainty 

Description 

Interval 

Gas 

He-40Ar ± 2% 

Gas composition 

0.02 

E 

2.0331 ± 6.1 xIO -4 V 

Terminal voltage 

0 

1 

569 ±0.74 A 

Current through test section 

0.001 

R 

3.58x10- 3 ± 4.77x10-6 ohm 

Resistance of test section 

0.001 

q 

1,156 ± 4.85 W 

Power supplied to inner tube 

0.004 

% 

-4.6 ±0.81 W 

Heat loss (insulation only) 

0.175 

<f" 

2.75x 10 7 ± 9.35x 10 5 W/m 3 

Volumetric heating rate 

0.034 

T t 

549.04 ±1.1 K 

Embedded thermocouple temperature 

0.004 

T w 

546.62 ± 1.95 K 

Surface temperature at tube OD 

0.008 

T 0 

311 ±1.1 K 

Outlet temperature from test section 

0.004 

T a 

300 ±1.1 K 

Room temperature 

0.004 

T 9 

291 ± 1.12 K 

Bulk gas temperature 

0.004 

T, 

271 ±1.1 K 

Inlet temperature to test section 

0.004 

^w'Tg 

255.62 ± 2.24 K 

Wall to gas temperature difference 

0.009 

TfT w 

2.43 ± 1.6 K 

Thermocouple to surface drop 

0.662 

L, 

25 ±0.01 K 

Test section total length 

0 

L a 

18 ±0.01 K 

Test section active length 

0.001 

X 

9 ± 0.1 in 

Measurement position 

0.011 

D s 

6 ±0.1 in 

Insulation outside diameter 

0.017 

D 0 

1.75 ±0.005 in 

Jacket outside diameter 

0.003 

Di 

1.25 ±0.002 in 

Jacket inside diameter 

0.002 

D, 

1.25 ±0.005 in 

Jacket total indicated runout 

0.004 

d o 

0.625 ± 0.005 in 

Inner tube outside diameter 

0.008 

d o 

0.625 ± 0.009 in 

Inner tube total indicated runout 

0.014 

d 

0.555 ± 0.005 in 

Thermocouple diameteral location 

0.009 

d i 

0.457 ±0.015 in 

Inner tube inside diameter 

0.033 

d h 

0.625 ± 0.004 in 

Hydraulic diameter 

0.006 

R 

0.084 ±0.016 in 

Inner tube wall thickness 

0.188 

r 

0.049 ±0.016 in 

Thermocouple to adiabatic face distance 

0.323 

A 

35.33 ± 0.344 in 2 

Heat transfer area 

0.01 

k 9 

0.066 ±2.5x10- 3 W/m-K 

Gas thermal conductivity 

0.038 

C 9 

1,131 ±1,13x10 1 J/kg-K 

Gas specific heat 

0.01 


2.28x10-5 ± 2.85 x10- 7 kg/m-s 

Gas viscosity 

0.013 

c s 

400 ±4x10 J/kg-K 

Stainless steel heat capacity 

0.01 

k s 

17 ± 1.7x1 0 _1 W/m-K 

Stainless steel thermal conductivity 

0.01 

Ps 

7,900 ±7.9x10! kg/m 3 

Stainless steel density 

0.01 

r s 

7.2 x 10“ 7 ±7.2x1 0- 9 ohm-m 

Stainless steel resistivity 

0.01 

k i 

0.068 ±6.8x10- 4 W/m-K 

Insulation conductivity vermiculite 

0.01 

e 

0 ± 0.032 

Eccentricity 

0.016 

Re 

30,000 ±1,294 

Reynolds number 

0.043 

Pr 

0.39 ± 0.02 

Prandtl number 

0.04 

Nu 

47.8 ±2.09 

Nusselt number 

0.044 
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